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ABSTRACT 

o ' 

CO ' We present the first direct measurement of the galaxy-matter bispectrum as a 

function of galaxy luminosity, stellar mass and SED type. Our analysis uses a galaxy- 
galaxy-galaxy lensing technique (G3L), on angular scales between 9 arcsec to 50 ar- 
cmin, to quantify (i) the excess surface mass density around galaxy pairs (excess mass 
, hereafter) and (ii) the excess shear-shear correlations around single galaxies, both of 

5^ ■ which yield a measure of two types of galaxy-matter bispectra. We apply our method 

to the state-of-the-art Canada- France-Hawaii Telescope Lensing Survey (CFHTLenS), 
spanning 154 square degrees. This survey allows us to detect a significant change of the 
bispectra with lens properties (stellar mass, luminosity and SED type). Measurements 
for lens populations with distinct rcdshift distributions become comparable by a newly 
devised normalisation technique. That will also aid future comparisons to other surveys 
or simulations. A significant dependence of the normalised G3L statistics on luminos- 
ity within -23 M r —18 and stellar mass within 5 x 10 9 M Q ^ Af* ^ 2 x 10 n M Q 
is found (h = 0.73). Both bispectra exhibit a stronger signal for more luminous lenses 
or those with higher stellar mass (up to a factor 2-3). This is accompanied by a steeper 
equilateral bispectrum for more luminous or higher stellar mass lenses for the excess 
mass. Importantly, we find the excess mass to be very sensitive to galaxy type as 
recently predicted with semi-analytic galaxy models: luminous (M r < —21) late-type 
galaxies show no detectable signal, while all excess mass detected for luminous galaxies 
seems to be associated with early-type galaxies. We also present the first observational 
constraints on third-order stochastic galaxy biasing parameters. 
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1 INTRODUCTION 

Over the course of the last two decades, the gravitational 
lensing effect has allowed us to establish a new branch of sci- 
ence that exploits the distortion of light bundles from distant 
galaxies ( "sources" ) in order to probe the large-scale gravi- 
tational field produced by intervening matter. Strong tidal 
gravitational fields cause an obvious distortion of individ- 
ual galaxy images ("strong lensing"; cf. Meylan et al. 2006), 
whereas weak deflections can only be inferred by statisti- 
cal methods utilising many galaxy images ("weak lensing"; 
cf. Schneider 2006). For the latter, usually shear image dis- 
tortions are harnessed, although the study of higher-order 
flexion distortions may also be feasible in the near future 
(cf. Goldberg & Natarajan 2002; Goldberg & Bacon 2005; 
Velander et al. 20f f ). Recently, the lensing magnification ef- 
fect has also moved into the focus of research as new source 
of information on cosmological large-scale structure (Hilde- 
brandt et al. 2009). As the gravitational field is solely de- 
termined by the mass density of the objects under exami- 
nation, no further assumptions on their properties need to 
be made when studying lensing. This makes it a unique tool 
for cosmologists to examine the large-scale structure of the 
Universe, in particular the relation between luminous com- 
ponents, such as galaxies, and the dark component. Within 
the current ACDM standard model of cosmology (Peacock 
f999; Dodelson 2003), the major fraction of matter is so- 
called dark matter, whereas ordinary baryonic matter is sub- 
dominant (Komatsu et al. 20fl). Therefore, lensing plays a 
key role in scrutinising the dominant matter component or 
in testing the standard model. 

Statistical methods have been developed that quan- 
tify the average mass distribution around galaxies by cross- 
correlating tangential shear, as observed from background 
sources, with foreground lens galaxy positions. Galaxy- 
galaxy lensing (GGL), as the first highly successful applica- 
tion, in effect measures the stacked projected surface mass 
density profiles around galaxies (Brainerd et al. f996; Hud- 
son et al. f998; Fischer et al. f999; McKay et al. 200f ; Hoek- 
stra et al. 2003; Sheldon et al. 2004; Seljak & Warren 2004; 
Hoekstra et al. 2004; Kleinheinrich et al. 2006; Mandelbaum 
et al. 2006; Parker et al. 2007; van Uitert et al. 201f ; Man- 
delbaum et al. 20f2; Leauthaud et al. 2012). The GGL signal 
is thus a function of lens-source separation (and their red- 
shifts) only, i.e., a two-point statistic that is based on a lens 
and the image ellipticity of a source galaxy. For a review see 
Schneider (2006) or Hoekstra & Jain (2008). GGL studies 
revealed, e.g., a mass distribution far exceeding the exten- 
sion of visible light: lenses are embedded in a dark matter 
halo of a size with at least ~ 100/i _1 kpc (Hoekstra et al. 

2004) and a mean density profile consistent with those found 
in ACDM simulations (Navarro et al. 1996; Springel et al. 

2005) . As extension of GGL, the light distribution within 
the lens can be utilised to align the stacked mass fields, 
which allows the measurement of the mean ellipticity of the 
halo mass distribution in a coordinate frame aligned with 
the stellar light distribution of a lens (Hoekstra et al. 2004; 
Mandelbaum et al. 2006; van Uitert et al. 2012; Schrabback 
& CFHTLenS team 2012). More generally, on larger spatial 
scales the technique has been exploited to infer the spatial 
distribution of lenses with respect to the matter distribu- 
tion, the second-order galaxy biasing (Hoekstra et al. 2001; 



Hoekstra et al. 2002; Pen et al. 2003; Sheldon et al. 2004; 
Seljak et al. 2005; Simon et al. 2007; Julio et al. 2012). More 
recently, GGL in combination with galaxy clustering in red- 
shift surveys has been employed to test general relativity 
(Reyes et al. 2010), or to successfully constrain cosmological 
parameters (Mandelbaum et al. 2012). 

[SW05] Schneider & Watts (2005, SW05 hereafter) in- 
troduced two new GGL correlation functions that involve 
three instead of two galaxies, either two lenses and one 
source ( "lens- lens-shear" ) or two sources and one lens ( "lens- 
shear-shear" ) . Therefore, this new class of correlators repre- 
sents the third-order level of GGL or simply "G3L". Both 
correlators express new aspects of the average matter dis- 
tribution around lenses, which can be translated into third- 
order galaxy biasing parameters (S W05) , especially if repre- 
sented in terms of aperture statistics (Schneider 1998). This 
paper chooses the aperture statistics to represent the G3L 
signal. Thereby we essentially express the angular bispec- 
trum of the (projected) matter-galaxy three-point correla- 
tion. A rigorous mathematical description of the aperture 
statistics is given in the following section. 

A more intuitive interpretation (Simon et al. 2012) of 
G3L is given by the definition of the real-space correlation 
functions: the lens-lens-shear correlation function measures 
the average excess shear (or excess mass, Simon et al. 2008) 
around clustered lens pairs, i.e., in excess of the average 
shear pattern around pairs formed from a hypothetical set 
of lenses that is uniformly randomly distributed on the sky 
(unclustered) but exhibit the same GGL signal as the lenses 
in the data. It is a probe for the joint matter environment of 
galaxy pairs, not single galaxies. This correlator promises to 
put additional constraints on galaxy models (Saghiha et al. 
2012) as it appears to be very sensitive to galaxy types. 
On the other hand, the lens-shear-shear correlation function 
measures the "excess shear-shear correlation": it quantifies 
the shear-shear correlation function in the neighbourhood 
of a lens in excess of shear-shear correlations as expected 
from randomly scattered lenses. Thereby it picks up the 
(projected) matter density two-point correlation function 
of matter physically associated with lenses. In a way this 
makes the lens-shear-shear correlator similar to the tradi- 
tional GGL, but now also probing the variance in the surface 
matter density around lenses instead of merely the average. 
The angular matter-galaxy bispectra are Fourier-transforms 
of these correlators. 

Simon et al. (2008) have demonstrated with the Red- 
Sequence Cluster Survey (RCS1; Gladders & Yee 2005) that 
both G3L correlation functions can readily be measured with 
existing lensing surveys. The RCS1 study aimed to obtain 
a high signal-to-noise ratio of the lensing signal, for which 
all available lenses were combined into one lens catalogue. 
Therefore, apart from this feasibility study in existing data, 
little more is known on the dependence of the G3L signal on 
galaxy properties. This paper is a first step to fill this gap 
by systematically measuring G3L for a series of lens sam- 
ples with varying properties. The amount of data available 
through the CFHTLenS analysis allows this to be done for 
the first time. An accompanying paper by Velander et al. 
(2012) explores the GGL signal of CFHTLenS in the light 
of the halo model (Cooray & Sheth 2002). 

The paper is laid out as follows. Sect. 2 summarises the 
aperture statistics that is devised to express the G3L signal, 
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gives their practical estimators and lists possible sources of 
systematics. In Sect. 3, we outline the selection criteria of 
our source and lens samples. Lenses are selected by luminos- 
ity, stellar mass, redshift and two galaxy spectral types, all 
to be analysed separately. Sect. 4 presents our G3L results. 
For a large range of angular scales covered in this study, the 
G3L signal is characterised by a simple power law whose pa- 
rameters are given. Sect. 5 offers a physical interpretation of 
the G3L statistics in terms of 3D galaxy-matter bispectra. 
In this context, we also introduce a normalisation scheme 
to remove, to lowest order, the impact of the exact shape of 
the lens redshift distribution and the source redshift distri- 
bution from the signal. Finally, the Sects. 6 and 7 present 
our discussion and conclusions. 

Throughout the paper we adopt a WMAP7 (Komatsu 
et al. 2011) fiducial cosmology for the matter density f2 m = 
0.27, the cosmological constant Qa = 1 — = 0.73 (both in 
units of the critical density) and Ho = 100/ikms~ 1 Mpc _1 . 
These parameters are consistent with gravitational lensing 
constraints obtained from CFHTLenS itself (Kilbinger et al. 
2012; Benjamin et al. 2012; Heymans et al. 2012). If not 
stated otherwise, we explicitly use h = 0.73, in particular 
for the absolute galaxy magnitudes and their stellar masses. 



2 FORMALISM 

This section summarises the theory and notation of G3L as 
detailed in SW05, and lists possible G3L specific systemat- 
ics. 

2.1 Galaxy-galaxy lensing preliminaries 

The weak gravitational lensing effect (see Schneider 2006, 
and references therein) probes the three-dimensional rela- 
tive matter density fluctuations 5 nl (R±,x) = A/9 m /p m in 
projection along the line-of-sight in terms of the lensing con- 
vergence 

Here R± = /k(x)^ is a 2D vector perpendicular to a refer- 
ence line-of-sight and the angular position on the sky. The 
comoving angular diameter distance /k(x) i s written as a 
function of comoving radial distance \- By Dh := c/Ho we 
define the Hubble length, and a(x) is the cosmic scale factor 
at a distance \\ we set a(0) = 1 by definition; c is the vacuum 
speed of light. By Xh we denote the comoving Hubble radius 
of today as the theoretical maximum distance at which we 
can observe objects. The lensing efficiency averaged over the 
probability density distribution function (p.d.f.) f>b(x)dx °f 
background galaxies ("sources") is expressed by 

9 (x)= r*x'Mx')^^- (2) 

Although the convergence in principle is observable 
through magnification of galaxy images, past weak lensing 
analyses and this paper focus on the related gravitational 
shear (Kaiser & Squires 1993) 

7c(0) = IJ - 0)k(4) ; V(0) := . (3) 



By 0* we denote the complex conjugate of 0. For this pur- 
pose, the complex ellipticity of the galaxy image 

e(0) ^ -y c (0) + e s ; (e s ) = (4) 

serves as a noisy estimator of 7 C ; the noise term originates 
from the unknown intrinsic shape e s . In addition, due to the 
finite number of sources, one also experiences sampling noise 
of the shear field. Note that we adopt the commonly used 
complex notation of 2D vectors and spinors (in the case of 
shears and ellipticities), where real and imaginary parts are 
the components along two Cartesian axes in a tangential 
plane on the sky. 

Galaxy-galaxy lensing techniques correlate the total 
matter distribution k(0) with the relative number density 
distribution k s (0) of lens galaxies ("lenses") on the sky by 
means of cross-correlating the lensing signal with positions 
of foreground galaxies, 

{e) = n g (0)-"E = T h d ^ sJficMO, x) , (5) 
«g Jo v ' 

where Pf(x)dx i s the p.d.f. of the lens (foreground) comov- 
ing distances along the line-of-sight; n s (0) is the projected 
number density of lenses and n s its statistical mean. For the 
scope of this paper, Pf(x) is estimated from a redshift p.d.f. 
p z (z)dz = p{ (x)dx of a selected lens sample. 

2.2 G3L aperture statistics 

For practical purposes, the aperture statistics are a conve- 
nient measure for a lensing analysis (Schneider et al. 1998; 
Schneider 1998; van Waerbeke 1998; Crittenden et al. 2002). 
They quantify moments of fluctuations in k(8) and n s (0) 
within apertures of a variable angular scale 9. The moments 
are determined from the smoothed fields k(0) and k s (0), 

M ap (0) = j^-u^O- 1 ) , (6) 

W = / ^-u^e- 1 ) (7) 

where u(-d/6)6~ 2 is the smoothing kernel. For mathematical 
convenience, we placed the aperture centre at C = in 
the previous definition. Third-order moments are defined by 
considering the ensemble average of 

(AA 2 M ap )(6»i;6»2;6'3) := <V(0i)Af(0 2 )M ap (0 3 )^ , (8) 

(NMl p ){6r,e2; 6 3 ) := /Af(ft)M ap (02)M ap (0 3 y , (9) 

over all random realisations of the fields k(0) and « g (0). Due 
to the assumed statistical homogeneity of the fields, the av- 
erages do not depend on the aperture centre position. There- 
fore, in practice, where only one realisation or survey is avail- 
able, these quantities are estimated by averaging the prod- 
ucts AA(6>i)Ar(6>2)Af ap (e 3 ) and AA(6»i)M ap (6l2)M a p(6l3) for dif- 
ferent aperture centres covering the survey area. See Fig. 1 
for an illustration. 

For a compensated filter u, i.e., J °° A8 8u{9) = 0, the 
aperture mass can in principle be obtained directly from the 
observable shear through (Schneider et al. 1998) 

M ap (0) = J™ fj q (M- 1 ) 5ft( 7 (tf ; v)) , (10) 
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Figure 1. Illustration of the aperture statistics. Fluctuations 
M(8) in the projected galaxy number density (top panel), 
smoothed to the characteristic filter scale 8, are statistically 
compared to the filtered projected matter fluctuations M ap (0) 
(lensing convergence; bottom panel). We take A^ 2 (#)M ap (0) or 
M (8) M%L(ff) , and average these for different aperture centres 
(dashed circles) to estimate third-order moments of the joint 
probability distribution of Af(8) and Af ap (#). 

where 7(1?;^) := — e~ 2l¥> 7 c (?9) denotes the Cartesian shear 
7 C at angular position i9 rotated by the polar angle <p. The 
real part of 7(1?; ip) is the tangential shear, the imaginary 
part the cross shear. The relation between the niters u(x) 
and q(x) is given by 

q{x)=(^r;j^ dssu(s)j -u(l) . (11) 

This paper uses the exponential aperture filter from van 
Waerbeke (1998), exponential filter hereafter, 

^)=^(i-£) e ~* 2/2 < < i2 > 

which effectively has a finite support because of the Gaussian 
factor that suppresses the filter strongly to zero for 1? > 30 
(SW05). The Fourier transform of the aperture filter is 

u(i) = J d 2 6 u{6)e +il - e = ^e" f2/2 . (13) 



Figure 2. Illustration of the paramctrisation of the lens- 
lcns-shear three-point correlator 6(t?i, 1?2 5 <p3) {t°P panel), and 
the lens-shear-shear correlation function G±(i?i , i?2i <p3) (bottom 
panel). These statistics are employed to estimate the aperture 
statistics in Fig. 1. The figure is copied from SW05. 

We generally denote a Fourier transform of f(0) by f(€) in 
the following. The exponential filter u(£) peaks in Fourier 
space at an angular wave number of £ = y/2, which deter- 
mines a characteristic angular scale selected by an aperture 
radius of 8. 



2.3 Aperture statistics estimators 

To obtain the third-order moments of the galaxy-matter 
aperture statistics, we utilise the lens-lens-shear correlation 
function Q in the case of (A/" 2 M ap ) and the lens-shear-shear 
correlation function G± for (A/"M ap ). This section provides 
only a brief description of this approach. For a more de- 
tails, its computationally optimised implementation as well 
as verification, we refer the reader to Sect. 3 of Simon et al. 
(2008). 

In practice, the aperture moments (A/" 2 M ap ) or {AfM^ p ) 
are not computed from the aperture mass M ap or aperture 
number counts N directly. The information contained in 
the aperture statistics is also contained inside two classes 
of three-point correlation functions (SW05), which are rel- 
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atively straightforward to estimate. Once the correlation 
functions have been determined, they can be transformed to 
the corresponding aperture statistics by an integral transfor- 
mation. The estimation process thus proceeds in two basic 
steps. In the first step, for (A/" 2 M ap ) one estimates the source 
tangential ellipticity relative to the midpoint connecting two 
lenses. 



e(tfi,tf 2 ,</>3) = i(n g (0i)n g (0 2 )Y(03; 



(14) 



The meaning of the notation is illustrated in the left panel 
of Fig. 2. For (NM 2 p ) one estimates the correlation of the 
cllipticities of two sources relative to the line connecting 
the sources as a function of separation from one lens (right 
panel), 

G±(0i,0a,6O = i-/7(0i;¥'i)7 ± (02;V2)%(03)\ • (15) 



Here and in the following equations a superscript "±" as 
in 7 ± means 7 for 7" (in case of G_) and the complex 
conjugate 7* for 7+ (in case of G+). 

Both correlation functions are estimated inside bins of 
similar triangles, i.e., lens-source triples within a configura- 
tion of comparable side lengths #1,2 and opening angles 03, 
by summing over all relevant galaxy triplets. Any triple of 
three galaxy positions 0i,0j,0k that meets the criteria of a 
relevant triangle is flagged by A^"' 2 ''' 3 = 1 and A^ 1 ^ 2 ^ 3 = 
otherwise. For this study, we utilise 100 logarithmic bins for 
both i?i and #2, and 100 linear bins for the opening angle 
</>3. For estimating Q we utilise 



N d N d N, 



(16) 



Jv d Jv d jv s r 1 

E E E Me- 1 ^^) 1 + 0,(10,-^1) Af. 

i=lj=lfc=l L 
JVd N d N B 

E E E 

i=l j=l fc=l 



k 



and for G± the estimator 



(17) 



JVd JVa N. 

E E I 

i=l j=l k=l 



E E E ™ 3 ™ fe e,e±e- 2i ^e ±2i ^A^ 2 ^ 3 



ijk 



E E E w w k At M 



ijk 



N d JV S JV B 

E E E 

i=l j=l k=l 

where Nd and N s are the number of lenses and sources, Wi 
are statistical weights of sources, ipi are polar angles of the 
position vectors of galaxies with respect to the coordinate 
origin, £j are the source ellipticities, and 



U (\A0\) = ( ^(0)^(0 + AO] 



(18) 



is the angular two-point clustering of the lenses (e.g. Pee- 
bles 1980). In this paper, the angular clustering of lenses 
is estimated by means of the estimator in Landy & Szalay 
(1993) prior to the estimation of Q and then interpolated. 
Sources are weighed by the inverse-variance uncertainty in 
the Zensfit ellipticity measurement (Miller et al. 2012). 

In a second step, we transform the estimates of Q and 
G± to the aperture statistics by devising the transforma- 
tion integrals Eqs. (63), (57), and (59) in SW05. There is 
no need to remove the unconnected terms in the correlation 



functions. As shown in SW05 (Sect. 7.2. therein), the trans- 
formation from Q to (A/" 2 M ap ) yields the same result when 
Q is taken instead of Q. Therefore, the integral transforma- 
tion automatically ignores unconnected second-order terms 
in the triple correlator, resulting in an aperture statistics 
that are only determined by pure (connected) third-order 
correlation terms. The same holds true for G± and (ATM™). 



2.4 Relation to 3D galaxy- matter bispectra 

The aperture statistics are directly connected to the angular 
cross-bispectra of the projected matter and lens distribution: 



{N 2 M ap )(6 1 ;6 2 ;8 3 ) 
d 2 h f d 2 l 2 



(2tt) 2 J (2tt) 



(UM! p )(e 1 -,e 2 -e 3 ) 



(19) 

■U(tiei)u(i202)u(\£l +l2\0 3 )bgg K (£l,l2) , 



(20) 



d 2 £i f d 2 £ 2 . 



)a u(*0i)ti(i2fc)ti(|*i +la\es)bn HS (li,l a ) , 
where the angular galaxy-galaxy-matter bispectrum is 
<j%(*i)Sg(*2)/e(40) = (27r) 2 ^ 2) (^+^+£ 3 )6g gK (^i^2)(21) 
and the angular matter-matter-galaxy bispectrum is 
(R(*i)is(Za)«g(40) = (27T) 2 5^ ) (e 1 +e 2 +£s)b KKg (e 1 ,£ 2 ) .(22) 

For statistically homogeneous random fields, the triple cor- 
relators on the left-hand side of the previous two equations 
can only be non-vanishing when t\ + t 2 + Is, — 0, which 
is reflected by the 2D Dirac delta functions 5^ (a;) on the 
right-hand sides. Owing to homogeneity, the bispectra thus 
depend only on two independent arguments £, for which 
we arbitrarily choose t\ and £ 2 . This automatically implies 
£3 = —(£1 + £ 2 ). In addition the statistical isotropy implies 
that the bispectra are solely functions of the moduli of £i i2 
and the angle enclosed by both wave vectors. 

As can be seen from Eqs. (19), (20), the aperture statis- 
tics are a locally filtered version of the bispectrum because 
the exponential u-filter is relatively localised in ^-space with 
a filter maximum at ^ ma x = v2/0. By means of the filtering, 
the aperture statistics basically becomes a band power bis- 
pectrum version of t> ggre or 6 KKg . Hence the aperture statistics 
Eqs. (8), (9) measure two different angular galaxy- matter 
band power cross-bispectra. 

By virtue of the Limber approximation (Kaiser 1992; 
Bartelmann & Schneider 2001) the angular bispectra and 
thereby the aperture statistics Eqs. (19), (20) can directly 
be related to the 3D cross-bispectrum of the matter and lens 
distribution (SW05) as primary physical quantities that are 
assessed by the statistics: 



b SSK (£i,e2) = 



3^ m r Xb , g(x)Pf (x) o ( £1 £2 \ 
• ' 3 ^^ ggm U(x)\/ K (x)' x J 



2£> H Jo fUxHx) 

b KKg (£i,£ 2 ) = 

9Ql fx* g 2 (x)Mx) p 



(23) 



(24) 



£1 



/k(x)' /k(x)' 



where the 3D bispectra are determined by the Fourier trans- 
forms of the matter density contrast, S m (k,x), and galaxy 
number density contrast, <5 g (fc,x), at radial distance x> 
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namely 

(d g (k 1 ,x)S g (k 2 ,x)S m (k 3 ,x))= (25) 

(2ir) 3 6g\ki + fc 2 + fc 3 )Bgg m (fc 1 ,fc 2 ,x) ' 
(Mfci,xMm(fe2,x)4(fe3,x)) = (26) 
(27r) 3 ^ 3) (fc 1 +fc 2 + fc 3 )B mmg (fei,fc2,x) ■ 

The vector k is the comoving wave number of modes entering 
the triple correlator. As before with the angular bispectra, 
the spatial bispectra are also isotropic, i.e., they are only 
functions of |fci|, &2| and the angle spanned by fci and &2- 
To refine the previous RCS1 measurement in Simon 
et al. (2008) for different galaxy populations, we focus 
on equally-sized apertures with 9i = 62 = #3 only. 
This leads us to the short hand notations {J\f 2 A4 a , P } (6) := 
(A/" 2 M ap )(6l;6l;6>), likewise for {AfM 2 p }. Due to the action 
of the it-filter in the Eqs. (f9) and (20) this picks up 
mainly bispectrum contributions from equilateral triangles 
l^i I = 1-^2 1 = |^i+^2 1, albeit also mixing in signal from other 
triangles because of the finite width of the u-filter in ^-space. 

2.5 Systematics indicators 

The gravitational shear of distant galaxy images is produced 
by small fluctuations 8(f) in the intervening gravitational po- 
tential. To lowest order in 3<f)/c 2 this is expected to only 
produce curl- free shear fields (B-modes vanish) . Current sur- 
veys do not have the power to measure higher-order effects, 
such that we expect these to be undetectable in our data. 
Shear-related correlation functions, or aperture moments in- 
volving the aperture mass, hence vanish after rotation of all 
sources by 45°, i.e., after 7 C (#) >-> — i7 c (0). Translated into 
data analysis, a 45° rotation of the source ellipticities should 
result in a measurement that is statistically consistent with 
the experimental noise (e.g. Hetterscheidt et al. 2007). We 
use this as a necessary (but not sufficient) indicator for the 
absence of systematics in the data. 

The estimator G± l in Eq. (17) incorporates two sources 
with two uniquely different possibilities to probe systemat- 
ics: rotating the ellipticities and tk of both sources results 
in the so-called B-mode channel of {AfM 2 p }(8), denoted here 
by (NMl)(6), and the P-mode channel, (NM±M ap }(0), if 
only either ej or are rotated. As pointed out by Schneider 
(2003), a P-mode is a signature of a parity- invariance viola- 
tion in the shear data, which in a parity-invariant universe 
can only be generated by systematics in the PSF correction 
pipeline, or in the algorithm for the statistical analysis of the 
data. Non-vanishing B-modes, on the other hand, can have 
a physical cause. For example, they can be associated with 
the intrinsic clustering of sources (Schneider et al. 2002) , in- 
trinsic alignment correlations of physically close sources or 
intrinsic shape-shear correlations (Heymans et al. 2006, and 
references therein). Especially the latter two are a concern 
for this analysis, as these effects are known to affect the E- 
mode channel of the aperture statistics, which is the prime 
focus of this work. However, currently it is unclear by how 
much this really affects G3L. We discuss in the following 
Sect. 2.6 that the influence of these systematics can be sup- 
pressed by separating lenses and sources in redshift, which 
is carried out in our analysis. 

Since the estimator Q est in Eq. (16) involves one source, 
there is only a single systematics indicator of {M 2 M SbP ){6), 



which is a parity violation indicator, a P-mode channel. In 
the following we will denote these statistics as {JV 2 M±}(8). 
As shown in SW05, the B- and P-modes of the statistics can 
be computed from Q and G± directly by utilising an alter- 
native integral kernel in the transformation from correlation 
functions to aperture statistics; see their Sect. 7.1 and 7.2. 

2.6 Reduction of II- and Gl-contributions 

One possible source of systematics are correlations with in- 
trinsic ellipticities e s of sources. A correlation between e a 
of different sources (Il-correlations) or between e s and a 
fluctuation in the mass density field generating shear (GI- 
correlations) is known to contribute to the shear correlation 
functions (e.g Hirata & Seljak 2004; Heymans et al. 2006; 
Joachimi et al. 2011). For a discussion of intrinsic align- 
ments in CFHTLenS see also Heymans et al. (2012). We 
argue here that selecting lenses and sources from well sep- 
arated distances ideally removes contaminations by II- or 
Gl-correlations in the G3L statistics. 

Consider the galaxy number density contrasts K gi i and 
K gi 2 in two arbitrary line-of-sight directions 6\ and 62, re- 
spectively, and a source ellipticity e s + 7 in a third direction 
63. The shear 7 and e s are rotated in direction of the mid- 
point between the two lenses according to the definition of 
Q. If lenses and sources are well separated in distance, then 
their properties are statistically independent. The lens-lens- 
shear correlator measures 

Q = (Kg,iK g ,2(7 + es)) (27) 

= (Kg,lK g ,27) + (Kg,lKg, 2 )(e s ) 
= («g,lKg,27) ) 

free of any systematic contribution from the intrinsic shape 
e s , if e s is statistically independent of the lens number den- 
sity fluctuation K g , i.e., 

(Kg,iK g ,2£s) = (/%,iK g ,2)(es) , (28) 

vanishing due to (e s ) = 0. 

Now, consider a lens number density contrast K g in one 
direction and the ellipticities e Sj , + 7i of two source images 
i = 1, 2 in two other directions. The ellipticities are rotated 
in direction of line connecting the sources in accordance with 
the definition of G±. The triple correlator measures 

G± = Mlt + ^ti)(l2 + e s , 2 )) (29) 

= vVYi 72) + (K g e^i72) + vVYi e f>.2> + (Kge s ,ie s , 2 ) 

= (regTFTa) + («g^i7a) + (^11^,2} + (««)(e^ie B ,2) 

= («g7i 7a) + (« g e s ,i72) + («g7i <*,2) ■ 

The last term in the third line vanishes because K g in 
the foreground is independent of the intrinsic shape of the 
sources in the background and because of (« g ) = 0. The 
latter follows from the definition of density fluctuations K g . 

The last two terms in the last line are less clear. For 
example in 172), 72 could be correlated with both e s ,i 

(GI signal, if source 2 is behind source 1) and K g (GGL 
signal). However, on the level of accuracy of the Born ap- 
proximation that is used in Eq. (1), the shear 72 is linear in 
the matter density contrast <5 m up to the distance of source 
2. We can, therefore, split the contributions to 72 into three 
parts 72 = 7« + 7 e + 7r OS t , namely (i) in contributions from 
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Figure 3. Prediction of the third-order moment of the aperture 
mass (M^ p ) for sources at redshift z B = 0.4. 

matter within correlation length to the lens, y K , (ii) matter 
within correlation distance to source 1, j t , and (iii) the rest 
7rest, which is neither correlated with k s nor with e B ,i- In 
this case we find 

(K g e^72> = («g7 K )(es ± i)+(«g><e 8 ± i7 e >+<« g )(e 8 ± i><7rest) -(30) 

All three terms vanish owing to (« g ) = (ejf^) = 0. A similar 
rational shows that also (f% 7j t e s ,2) vanishes to lowest order, 
such that we expect to find in the weak lensing regime 

G± = <« g 7i t 72) • (31) 
2.7 Magnification of lenses 

Another conceivable systematic effect is through cosmic 
magnification (Narayan 1989; Bartelmann & Schneider 
2001) that is generated by matter density fluctuations in 
front of lenses. To lowest order, foreground matter density 
fluctuations with lensing convergence k< (Eq. 1) integrated 
to the lens distance modify the observed clustering of lenses 
on the sky above a certain flux limit fy im according to 

k' s = K s + Xk< + 0(k%) , (32) 

compared to the unmagnified lens number density « g . Here, 
we have A := 2(v—l) with n g (> /ii m ) oc being the mean 
number density of lenses with flux greater than /n m . Nor- 
mally v — 1 is of the order of unity (van Waerbeke 2010) 
or smaller. Likewise the shear distortion 7 = 7< + 7>, 
Eq. (3), into the same l.o.s. direction contains a contri- 
bution 7< related to «:<, and 7> that is the shear orig- 
inating from matter fluctuations beyond the foreground. 
This in combination produces as additional contribution to 
S = (« g ,iK g ,2(7< + 7>)) the term A 2 (k<,ik<,27<) and to 
G ± = <Kg(7<.i+7>,i)(7<,2+7>.2)) theterm A(k<7< : i7<,2). 

These terms are basically third-order cosmic shear cor- 
relations or, in terms of the aperture statistics, related to the 
(Ma P (0)) statistics (Schneider et al. 2005). Third-order shear 
correlations have been measured (Bernardeau et al. 2003; 
Pen et al. 2003; Jarvis et al. 2004; Semboloni et al. 2011), 
and (M,L(#)) has been found (Jarvis et al. 2004; Semboloni 
et al. 2011) to be of the order of < 10~ 7 for aperture scales 



Figure 4. Total number of lenses (red) and sources (blue) in 
the catalogue between 0.2 ^ £ p hoto < 1-2 and 17.5 ^ i' < 22.5 for 
lenses or 17.5 ^ i' < 24.7 for sources. The figures comprise all 
galaxies, complying with the selection cuts, contained within all 
172 pointings. For the G3L analysis, the lens sample is further 
subdivided in luminosity, stellar mass bins and photometric red- 
shift, while sources are rejected for a photo-z of Zphoto < 0.65. 

of 6 ~ l' and sources at z s ~ 1.0. As this includes contribu- 
tions from the entire integrated matter up to z s , whereas the 
G3L magnification effect only contributions from the matter 
integrated up to the lens redshifts Zd ~ 0.4, we consider this 
an empirical upper limit for the magnification effect. In Fig. 
3, we show a prediction of (M& p (8)} with sources at z s — 0.4 
for a WMAP7-like cosmology based on the theory described 
in Semboloni et al. (2011). This result implies that the im- 
pact of lens magnification on the G3L aperture statistics is 
smaller than < 10 -8 . 



3 DATA 

3.1 Object selection and photometric redshifts 

This work uses the full CFHTLenS data set, which origi- 
nates from the CFHTLS-Wide Survey. The CFHTLS-Wide 
imaged 171 MegaCam (mounted on the CFHT) pointings in 
the five broad-band filters it*, g', r' , i' , and z' . During the 
observation campaign of CFHTLS, the j'-band filter was re- 
placed by a new filter with a slightly different transmission 
curve. For some of the pointings only the updated i'-band 
filter magnitudes are available, which are treated as the old 
filter magnitudes in the analysis. For details, see Erben et al. 
(2012). 

CFHTLenS has an effective area (different pointings 
partly overlap) of about 154 square-degrees with high- 
quality photometric redshifts down to il w 24.7. The data 
set and the extraction of our photometric redshift catalogue 
are described in Hildebrandt et al. (2012). Our data pro- 
cessing techniques and recipes are described in Erben et al. 
(2009) and Erben et al. (2012). As primary selection cri- 
terion, we select sources brighter than i' < 24.7 and lenses 
brighter than i' < 22.5. This will be further subdivided in 
the following by using photometric redshifts (Fig. 4) and, in 
the case of lenses, M r rest frame magnitudes, stellar masses 
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Table 1. Selection criteria of lens samples and source sample for the G3L analysis applied to the samples in Fig. 4, following Mandelbaum 
et al. (2006) for the lenses. The luminosity bins (L), stellar mass bins (sm) and galaxy type bins (ETG: early- type galaxies; LTG: late- 
type galaxies) are again subdivided by 0.2 < 2 p j loto < 0.44 ("low-2") and 0.44 ^ z p hoto < 0.6 ("high-2"). Sources attributed no statistical 
weight w by Zensfit are not used in the source sample. The galaxy numbers are for all pointings of which the final analysis discards 
roughly 25%. Luminosities and stellar masses assume h = 0.73. (1) z: mean redshift, cr z : r.m.s. variance of p(z); (2) and (3): best-fit 
parameters of w(8) = A u {6 /X')~ * + IC within 0'.2 < 9 < 10'; (4): sample completeness; (5): mean r-band luminosity; (6): mean stellar 
mass in units of 10 10 Mq. 
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or SED information (details below). 43 pointings out of 171 
exhibit a significant PSF residual signal, according to the 
detailed tests in Sect. 4.2 of Heymans et al. (2012), and are 
therefore discarded for the analysis (~ 25% area); 129 point- 
ings are included in the analysis. This leaves a total effective 
survey area of ~ 120 deg 2 that is eventually used in the anal- 
ysis. Of this area an additional ~ 20% percent is lost due 
to masking. The analysis is performed on individual fields 
which allows us to use field-to-field variances of the mea- 
surements to estimate the covariance of measurement errors 
directly from the data. 



3.2 Lens samples 

To guarantee a high reliability of the photo-z estimates for 
the lenses, a magnitude cut of i! ^ 22.5 is applied. A detailed 
account and tests of the CFHTLenS photo-z pipeline can 
be found in Hildebrandt et al. (2012). Based on the galaxies 



endowed with photometric redshifts, three classes of lens 
samples are selected (Table 1): 

• A luminosity or L-sample class, which consists of six 
distinct rest- frame M r bins (SDSS r-filter; York et al. 2000), 
labelled LI to L6. The same formal luminosity bin limits as 
in Mandelbaum et al. (2006) or Velander et al. (2012) are 
applied, although we do not automatically expect equivalent 
completeness of the samples. To quantify the completeness, 
we introduce the f c parameter below. 

• A stellar mass or sm-sample class, which is also further 
subdivided using seven distinct stellar mass bins. Again, we 
are guided by Mandelbaum et al. (2006) for compiling this 
sample class. The sm class has sub-classes with labels sml- 
sm7. 

• A galaxy type class using the T_B parameter in BPZ 
(Benitez 2000), which provides the most likely galaxy SED 
for a given galaxy and its estimated photo-z; see Erben et al. 
(2012) for more details. T_B=2 as division line, we separate 
early- type galaxies ("ETG"), which have T_B< 2, from late 
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Figure 5. In contrast to Fig. 4, this figure shows the full BPZ redshift posterior of the various samples. Low-z lenses are selected from 
2 photo £ [0.2,0.44], high-z lenses from z p hoto £ [0.44,0.6], and sources from z p hoto £ [0.65,1.2]. The dashed black line enclosing the 
source p.d.f. is a parametrised fit, see Sect. 3.3 for best-fit parameters. 



type galaxies ("LTG"). 1 In order to define a volume-limited 
sample of ETG and LTG, we select only luminous galaxies 
with restframe luminosities — 23 ^ M r < —21. With this 
luminosity cut, ETG and LTG are actually subsamples of 
L5 and L6 combined. 

The stellar masses of the lenses are determined from the 
galaxy multi-colour data as described in Sect. 2.1 of Velander 



Within BPZ values of T_B denote best-fitting galaxy tem- 
plates: 1=CWW-E11, 2=CWW-Sbc, 3=CWW-Scd, 4=CWW-Im, 
5=KIN-SB3, and 6=KIN-SB2. Note that the templates are inter- 
polated, such that fractional numbers occur. 



et al. (2012). The estimators assume a Chabrier (2003) star 
initial mass function. 



All three classes are further split into two photo- z bins: 
a "low-z" bin with 0.2 ^ 2 p hoto < 0.44 and a "high-z" bin 
with 0.44 ^ Zphoto < 0.60. As redshift estimators we use the 
maximum probability redshifts of the redshift posterior pro- 
vided by BPZ. The redshift boundaries give comparable num- 
bers of lenses prior to attributing them to one of the three 
lens classes (Fig. 4). Not counting the high-z LI and L2 sam- 
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pies, which have too faint limits to contain lenses 2 , we have 
in total 28 lens subsamples. 

The true redshift distribution of a lens sample is 
not identical to the distribution of their photometric red- 
shifts due to the errors in the photo-z estimators. For a 
magnitude cut of i' < 22.5, the errors are approximately 
°z ;$ 0.04(1 + z) with a ~ 3% outlier rate (Hildebrandt et al. 
2012). We combine the posterior redshift probability distri- 
bution functions (p.d.f.) of all lenses given by BPZ, see Fig. 5, 
to quantify the redshift uncertainties of complete lens sam- 
ples. The depicted redshift probability distributions will be 
utilised when normalising the G3L aperture statistics. 

To help the comparison of our G3L results to future 
studies, we also quote the angular clustering and complete- 
ness of the lens samples. The results are listed in Table 1, 
the details are described in Appendices Bl (clustering) and 
B2 (completeness). In short, for the angular clustering of 
lenses, we approximate the angular galaxy two-point cor- 
relation function by a power law over the angular range 
0'.2 ^ 9 < 10'. For each lens sample with the photo-z bin 
[zi , Z2] , we quote the completeness factor / c that expresses 
the average V(zx, z max )/V r (zi, Z2) of all lenses in the sample; 
V{z\,Z2) is the light cone volume between redshift z\ and 
Z2, and z max ^ Z2 is the maximum redshift up to which a 
lens is still above the flux limit i' = 22.5. A small f c is a sign 
of a strong incompleteness because many galaxies similar to 
those observed near z\ are presumably missing at higher 
redshifts. Due to the magnitude limit, samples containing 
a substantial portion of faint galaxies are most affected by 
incompleteness, most notably LI and L2. As expected, the 
completeness drops if one moves from the low-z to the high- 
z bin in almost all cases. The few minor exceptions, L5 for 
instance, are probably due to shot noise in the f c estima- 
tor. We conclude that L4-L6, sm3-sm7 and ETG/LTG are 
the most complete, volume-limited samples for our study 
(f c > 0.80 for both low-z and high-z). In Table 1 we also 
quote the average absolute r-band flux of the samples, listed 
as magnitude (M r ) and the their average stellar mass (M»). 

3.3 Source sample 

All details concerning the galaxy shape measurement (em- 
ploying the lensfit algorithm; Kitching et al. 2008, Miller 
et al. 2007, and Miller et al. 2012), CFHTLenS source cata- 
logue generation, and the discussion of shear systematics are 
presented in Heymans et al. (2012) and Miller et al. (2012). 
We account for the multiplicative shear bias by employing 
the Miller et al. (2012) normalisation scheme adjusted to 
our estimators (see Appendix A). 

In order to reduce the level of undesired II- and GI- 
correlations in the measurements, we attempt to separate 
sources and lenses by redshift, utilising photometric red- 
shifts as estimators. As a compromise between accurate red- 
shift estimates and a large numbers of sources, we apply a 
magnitude limit i' < 24.7 to the lensfit shear catalogue and 
select sources between 0.65 ^ z p hoto < 1-2. As for the lenses, 
the true redshift distribution is derived from the combined 

2 Actually, we find a few galaxies in the high-2 L1/L2 samples. 
These are probably extreme outliers with greatly inaccurate red- 
shift estimates. 



posterior redshift p.d.f. of individual sources, shown in ev- 
ery panel of Fig. 5 in comparison to the redshift distribution 
of the lens samples. The individual posteriors are weighted 
with the source weight that is also used in the lensing anal- 
ysis. The source redshift p.d.f. is well fitted by a broken 
exponential distribution 

PbW cJ exp ;- p f°-^;j if *^°< (33) 

[ exp (— p2(z — zo) Pa ) otherwise 

with fit parameters po = 91.14, pi = 2.623, P2 — 4.093, 
pz = 1.378 and z = 0.794 (dashed black lines). With our 
selections we find about 3 x 10 6 sources with mean redshift 
z ~ 0.93. As can be seen in Fig. 5, the overlap of the various 
Pf (z) and the source Pb(z) is small but not entirely vanishing, 
mainly at z = 0.5 — 0.7 for the high-z and at z ~ 0.6 for the 
low-z samples. The typical overlapping area of the redshift 
probability distribution functions (visible in Fig. 5) is ~ 12% 
for the high-z samples and ~ 4% for the low-z samples. 



4 RESULTS 

4.1 Measurements and their uncertainties 

In order to obtain measurements for the lensing aperture 
statistics, we use the method outlined in Sect. 2.3. As the 
binning grid for Q es and GJ 5 * , 100 log-bins ranging between 
9 arcsec and 50 arcmin are set up for #1/2, 100 linear bins are 
used for the opening angle 4>3, yielding overall 10 6 bins with 
bin widths A<^>3 = 3.6 deg and A In?? = 0.058. All measure- 
ments are performed separately on every individual point- 
ing, out of 129 square pointings with roughly 1 deg 2 each. 
Adjacent pointings partly overlap, however, which reduces 
the area that is actually used. In our study, we crop the 
pointings to remove the overlap. For the final result, indi- 
vidual estimates are combined by averaging the individual 
Q cst and G±* weighted by the number of triangles within 
each bin. 

Finally, the combined estimates are transformed to 
the aperture statistics by the integral transformations dis- 
cussed in SW05. In his way, the aperture statistics between 
0!5 ^ S ap ^ 10' for ten aperture scale radii are computed. As 
addressed in Simon et al. (2008), the transformation from Q 
or G± to aperture statistics becomes biased towards small 
and large aperture radii due to an insufficient sampling of 
the correlation functions. A similar transformation bias is 
also known for the aperture mass statistics (Kilbinger et al. 
2006). For the small separations, the bias depends in detail 
on the mean number density of the galaxies, most crucially 
the lenses, and the clustering of the lenses, which in com- 
bination determines the sampling of the correlation func- 
tions by small triangles. By comparison to simulated data, 
we made sure that this bias is negligible (below ~ 10%) 
within the range of 1' < 6 < 10' in our case (see Fig. 1 in 
Saghiha et al. 2012 for an illustration of the transforma- 
tion bias). The variance of the measurements across all 129 
pointings is used to estimate the covariance of measurement 
errors (Jackknifing; Appendix Bl). The inverse covariance 
matrix is estimated from the pointing-to-pointing covariance 
according to the method in Hartlap et al. (2007). 
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Figure 6. Measurements of the E-mode aperture statistics (M 2 M ap } (9) (top figure) and (AfM^ p )(9) (bottom figure) as a function of 
aperture scale radius 9. The left column depicts measurements for the low-z bin, the right column the high-z bin. Different lines refer 
to different lens samples (Table 1). Note that the values get biased for 9 < I' due to the transformation bias. Error bars indicate the 
la standard deviation of the mean of all pointings considered. Missing data points are outside the plotting range but consistent with 
zero. Numbers at the top indicate the effective scale of the statistics according to the maximum in the u-filter. L1-L6: M r -luminositics 
increasing from -17.8 mag to -22.4 mag; sml-sm6: increasing stellar masses from 7 X 10 9 Mq to 2 X 10 11 M . 
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Figure 7. Aperture statistics results for the late-type and early-type lens samples. The left column shows the low- ,2 rcdshift bin, the 
right column the high- z redshift bin. Top row and bottom row correspond to (A^ 2 M ap ) and (A^M 2 p ), respectively. Numbers at the top 
indicate the effective scale according to the maximum in the n-filter. ETG: early-type galaxies; LTG: late-type galaxies. 



Table 2. Results of the null tests for the different channels of the statistics. Quoted values are the reduced Ax 2 with d.o.f=7. for the 
E-modes and P-modc, {Af 2 Mj_), and d.o.f.=14 for the combined B/P-mode, (A/"Af?) and {A^MxM ap }. E-modes consistent with a null 
signal are in bold, P- or P/B-modcs inconsistent with a null signal arc underlined . Adopted confidence levels are 95%. Only data points 
within the range 1{5 ^ 9 < 10' were used for this test. 

0.22 sC 2photo < 0.44 0.44 z photo < 0.60 

<Af 2 A/ap> <AfM 2 p > <Af 2 A/ap> <AfM 2 p > 

Sample E-mode P-mode E-mode P/B-mode E-mode P-mode E-mode P/B-modc 



LI 


3.22 


1.27 


L2 


7.32 


1.57 


L3 


3.79 


0.66 


LI 


9.26 


0.58 


L5 


6.72 


1.00 


L6 


7.76 


0.99 


sml 


6.41 


1.06 


sm2 


12.43 


0.59 


sm3 


7.07 


0.99 


sm4 


7.64 


1.86 


sm5 


3.76 


1.27 


sm6 


0.65 


0.91 


sm7 


0.93 


1.16 


ETG 


12.50 


0.87 


LTG 


0.77 


0.85 



2.64 


1.24 




2.31 


1.75 




4.93 


0.75 


3.21 


2.62 


1.14 


7.66 


4.31 


0.80 


7.08 


4.92 


0.74 


7.07 


4.06 


1.55 


5.18 


3.55 


1.11 


5.97 


4.05 


0.76 


5.62 


4.83 


0.58 


5.93 


4.95 


0.52 


6.74 


1.15 


1.00 


2.33 


0.66 


1.61 


0.45 


6.62 


0.65 


15.24 


1.58 


1.52 


0.90 



0.94 


2.38 


1.22 


0.89 


2.38 


0.79 


1.95 


0.81 


1.06 


1.50 


1.74 


0.41 


0.74 


1.90 


0.68 


0.40 


0.95 


0.83 


0.65 


0.18 


1.19 


0.27 


2.28 


1.12 


1.18 


0.78 


0.96 


0.68 


0.64 


1.42 


0.42 


1.52 


1.87 


1.74 


1.21 


1.06 


1.20 


1.45 


0.37 
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4.2 E-mode measurements 

Fig. 6 summarises the E-mode results for the luminosity 
and stellar mass bins of the (A/" 2 M ap } (top) and (AfM 2 p ) 
statistics (bottom). Due to the incompleteness in the sam- 
ples, LI and L2 are empty in the higher redshift bin and 
hence are missing in the corresponding plots. Likewise, due 
to the small number of lenses and correspondingly large er- 
ror bars, also the data points of sm7 are missing. The signal 
dependence on galaxy type is displayed separately in Fig. 
7. For aperture radii greater than ~ 2 arcmin the measure- 
ments seem to be well approximated by power laws, which 
will be determined below. Below roughly 2 arcmin there are 
indications of deviations from the power-law behaviour at 
smaller radii in several cases, e.g., (AT 2 M ap } of low- z L1/L4, 
or (JVMfp) of high-z L4/L6. 

The result of (M 2 M a , p ) of the late-type galaxies (LTG) 
stand out as being the only one that is completely consistent 
with zero despite relatively small error bars. Therefore, the 
excess mass around late-type galaxy pairs vanishes within 
the statistical uncertainties. In strong contrast to that, the 
corresponding signal of the early-type galaxy (ETG) sample 
is highly significant. From the LTG signal upper limit we 
estimate the ETG signal to be greater by a factor of at least 
~ 10. This confirms the prediction of Saghiha et al. (2012) 
that is based on galaxy population synthesis models. 

The low- z sample LI, with the fewest number of lenses, 
presumably is affected by the transformation bias. This can 
be seen by the clear drop of the data points for (A/" 2 M ap ) 
below S ap ~ 2 arcmin compared to a power-law behaviour 
at larger scales. 

4.3 Systematics tests 

General tests for systematics on the level of shear catalogue 
generation are to be found in Heymans et al. (2012). We only 
use CFHTLenS pointings that passed the therein described 
tests for cosmic shear applications. To further test for sys- 
tematics in our measurements, we check for the consistency 
of the aperture statistics B- and P-modes with a null signal. 
The details of this test and, moreover, G3L measurements 
within separate CFHTLenS fields (W1-W4) are presented in 
Appendix C. The null test also allows to quantify the signif- 
icance of the signal in the E-mode channels of the statistics. 
Table 2 summarises the tests for all statistics and galaxy 
samples. 

In summary, we find that B/P-modes in the aperture 
statistics are consistent with zero between 1 arcmin and 10 
arcmin. When looking at the combined L1-L6 sample, sep- 
arate measurements within the survey fields W1-W4 agree 
well for both the low-z and the high-z redshift bin. This 
demonstrates the internal consistency of the data, and that 
the observed signals do not originate from a single, possibly 
peculiar field. As to the E-mode channels of the statistics, 
we find for (A/" 2 M ap } highly significant signals (95% confi- 
dence) for all low-z samples, except for sm6 and sm7, and 
most high-z lens samples. Sm6 and sm7 pose exceptions 
because they contain relatively small numbers of galaxies. 
Apart from the high-z L3, L4, and sm4, all high-z sample 
measurements of {JVM 2 p ) are consistent with zero, whereas 
their low-z counterparts are mostly significant. As {MM 2 p ) 
involves a three-point correlation function with two sources 
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Figure 8. Visualisation of the (^-measurement. Shown is the ex- 
cess mass (/cg(0i)ftg(02)t(03)) (intensity scale) as function of 0;j 
around early-type (top) and late-type galaxy pairs (bottom) with 
mean angular separations \0\ — 02 | between 40 and 60 arcsec; one 
map pixel corresponds to 1.67 arcsec (x- and y-axis labels). The 
lens positions 6\ , O2 are indicated as boxes, the contours show 
the S/N levels 2, 3, . . . (positive excess mass) and -3,-2 (negative 
excess mass). To increase the signal-to- noise the low-z and high- 
z maps of the ETG and LTG samples have been combined. A 
smoothing with a Gaussian kernel of r.m.s.-size 6.7 arcsec has 
been applied to map. 



and one lens, the noise level of this measurement is naturally 
higher than for (A/" 2 Af ap ). 
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5 INTERPRETATION 
5.1 Lens-pair excess mass 

Although aperture statistics and G3L-correlators essentially 
contain the same information, we would like to show our Q- 
measurements for at least the ETG and LTG samples. As 
outlined in Simon et al. (2008), the G3L correlation func- 
tion Q can conveniently be interpreted as a convergence map 
(excess mass map) once the separation of the two lenses is 
fixed; the E-mode in (J is a series of such maps for varying 
lens-lens separations. After a rotation, the correlator Q is 
a stacked shear field around the lens pair, from which we 
subtract off the GGL signal around individual the lenses to 
determine the connected part Q. To obtain the excess mass 
maps for the ETG and LTG in Fig. 8 we transform this 
stacked shear field to a convergence map utilising the algo- 
rithm in Kaiser & Squires (1993). For these maps, we con- 
sider relatively small lens-lens separations between 40 and 60 
arcsec as in Simon et al. (2008), and we combine the maps 
of the low-z and high-z samples; lens-lens-source triangles 
are rescaled inside the map such that lenses are always at 
the same position in the map (boxes). We also exploit the 
parity invariance of the maps by averaging the left and right 
half of the map, thereby increasing the signal-to-noise, see 
Simon et al. (2008) for details. 

The ETG map contains more significant structure and 
higher convergence values compared to the LTG map, which 
has only a weak signal. Qualitatively, the excess mass of the 
ETG sample is concentrated between the lens pair, whereas 
the LTG lenses seem to possess a small halo of excess mass 
around the individual lenses and a convergence trough be- 
tween them. The latter implies that the average convergence 
about a LTG pair (both lenses at similar distance) at given 
separation is lower than the sum of convergence around two 
mean individual late-type galaxies. We will study these maps 
in more detail in a forthcoming paper and focus on the aper- 
ture statistics for the remainder of this paper. 



5.2 Power-law fits to measurements 

For aperture scale radii larger than 9 > 2', our measure- 
ments are reasonably consistent with a simple power law. 
Therefore, we fit power laws (A/" n M™)(0 ap ) = A (6/l')- a 
to data points within l!5 ^ 9 < 10' to quantify the measured 
profiles of the statistics n = 2, m = 1 and n = 1, m = 2, see 
Table 3. The fit starts at 1.5 arcmin in order not to be too 
strongly influenced by the transformation bias. Fits use the 
Jackknife covariance matrices based on the measurements 
in the different pointings, as in Eq. (B2). For the fit, a mul- 
tivariate Gaussian noise model for the measurement errors 
is adopted. The quoted values indicate the posterior me- 
dian and a 68% credibility region about the median for the 
amplitude Ao and slope a. The posterior adopts a top-hat 
prior for the power-law slope, only allowing values within 
a € [0,5]. 

Fig. 9 depicts the dependence of the slope a on the lens 
M r magnitude and stellar mass for (AT 2 Af ap ) of all sam- 
ples with at least a 95% confidence detection (Table 2). We 
find a clear trend towards steeper slopes (steeper equilat- 
eral bispectra) for more luminous galaxies and galaxies with 
higher stellar mass. Note that sm6 and sm7 contain on av- 



erage galaxies more luminous than these of L5 and L6 (Ta- 
ble 1). The figure also depicts the measured slopes for the 
ETG samples, which are consistent with the L-subsamples 
of comparable M r luminosity (between L5 and L6). Slopes 
weakly constrained by the data have posterior medians that 
are drawn towards the centre of the top-hat prior, a — 2.5. 
This mainly applies to the noisier (MM 2 p ) measurements, 
for which reason they are not included in Fig. 9 but are 
listed in Table 3. 



5.3 Normalised measurements 

The G3L aperture statistics are directly related to the 3D 
matter-galaxy cross-bispectra and the redshift distribution 
of lenses and sources (Sect. 2.4). The radial galaxy distribu- 
tions and the fiducial cosmology define a smoothing kernel 
in radial and transverse direction. To disentangle, to low- 
est order, the dependence of the signal on the physically 
relevant bispectrum from the dependence on source or lens 
distribution, we introduce a normalisation scheme. 

Combining Eq. (19) and Eq. (23) shows that (A/" 2 M ap ) 
constitutes a radially and transversely weighted average of 
the 3D bispectrum _B ggm (fci, &2, x)i namely 



<ArM ap >(6>i;6>2;6>3) 

- ^ f XH A g(x)Pf(x) 

12 / Q A r 3 . 



(34) 
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By changing the integration variables as in £j = /kW^ we 
write this integral as 



<A/" 2 Af ap )(0 i; 02;03) 

dXfem(x) 



(35) 
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for which we introduce the tt-filtered bispectrum 
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As implied by (35), the lensing aperture statistics is ba- 
sically the transversely u-filtered bispectrum B g gm aver- 
aged in radial direction by the kernel Q ggm (x)- F° r equally- 
sized aperture radii 9\ =62 = #3, the ii-filter gives most 
weight to the equilateral bispectrum B ggm (fci, ft2, x) with 
fci = = |fei + fca|, but also mixes in other triangle config- 
urations to some degree. The radial weighting kernel, 



3n m 9(x)pUx)Mx) 

2 a( X ) 



(37) 



is peaked at a radial distance Xmax that is determined by 
the redshift p.d.f. of lenses and sources (top row of Fig. 10). 
Therefore, most weight is given to the bispectrum at dis- 
tance Xmax- 
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Figure 9. Dependence of the power-law index a in the normalised aperture statistics (A^ 2 M ap ) between 1'5 ^ 9 < 10' with M r -band 
magnitude (left panel) and stellar mass of lenses (right panel) . Table 3 lists the power law indices of all statistics. Only fits to measurements 
with a 95% confidence detection are plotted. The shaded region highlights the 68% credibility region of the combined low-z (open crosses) 
and high-z (open stars). Also shown are the slopes for the early- type galaxy sample in the corresponding magnitude range (left panel). 
For clarity, these data points are offset about their actual mean (M r ) = — 21. 88(— 21.91) for the low-z (high-z) sample. 



Table 3. Power-law fits (J\f n M™)(9) = Aq{8/1')~ a to the measured aperture statistics in Fig. 6; Aq is the signal amplitude at an 
aperture scale radius of = 1 arcmin. The fit considers only data within 9 £ [lf5, 10'] . Quoted errors bracket a 68% credibility region 
about the median. 

0.20 < z photo < 0.44 0.44 < z photo < 0.60 

<A/" 2 A/ap> <A/-A4 2 p> (A^ 2 Af ap ) {MMl p ) 

Sample Ao/lO" 5 a A /IQ- 7 a Ao/lO" 5 a A /1Q- 7 a 



LI 


10 


.-+5.12 
°' -4.90 


1.40 


L2 
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0.97 


L5 


12 


247 
30-2.33 


1.44 


L6 


35 


^4+9. 15 
-7.80 


1.91 



+0.90 
-0.34 
+0.18 
-0.16 
+0.31 
-0.24 
+0.20 
-0.17 
+0.23 
-0.18 
+0.35 
-0.26 



54.38 
9.48 



+25.92 
20.06 
+8.30 
11.37 



2.59"* 



no qq+7.95 
zo.zz_g.63 

iD '' D -6.88 
34 5K+14.41 

-+25.33 
18.24 



62.85] 



1+0.94 
-0.59 
n i 1+1-59 
z ' 11 -0.73 
1 C7+0.44 
l-»'_0.29 
1 95+ 1 - 14 

1 QO+0.80 
i - 9y -0.37 

2 in+°- 58 

z.1U_q 3 g 



8.56 
6.85 
9.15 
19.89 



+3.78 
2.75 
+ 1.50 
1.34 
1.80 
1.70 
+6.28 
5.48 



1.66 
1.49 
1.50 
1.68 



0.84 
0.41 
+0.29 
0.23 
+0.26 
0.20 
+0.31 
0.30 



^J-U' -21.95 
SI qo+ 18 06 

o oq+14.23 
-8.89_lg.56 

42 n7+ 29 - 43 



3.84 
3.30 
3.58 
3.43 



+0.71 
1.06 
+0.95 
0.87 
+ 0.88 
1.33 
+ 0.93 
1.02 



sml 
sm2 
sm3 
sm4 
sm5 
sm6 
sm7 



-+2.57 
-2.42 
,+3.09 

-2.96 
3 +5.22 
-4.91 
-,+6.88 
-6.40 

87.19133-76 
197.061313-36 
315 06+ 307 - 91 



11.95] 
19. 9lj 
24.49]" 
34.96"' 



1 OO+O- 22 
± - za -0.18 
4+0.17 
-0.14 
1 54+ - 18 
-+0.25 
-0.19 
i o 7 +0.65 
1 - 8 '-0.27 
3+0.84 
-1.02 
,+1.04 
1.32 



1.24 1 



1.65 1 



3.60] 
3.26^ 



42.70 
38.15 
45.06 
61.43 
77.89 
124.38 
230.78 



19.94 
12.91 
+ 16.87 
12.24 
l«.!ir. 
14.01 
+ 19.67 
16.56 
+26.49 
22.99 
+83.56 
64.46 
+265.46 
198.85 



2 16+ 095 
z ' iD -0.46 
3+0.82 
-0.39 

z,u3 -0.36 

1 95+ - 51 
l.»o_ 2g 

1 96+ 45 
J..»u_ 2g 

,+1.00 

-0.91 

,+0.74 

1.19 



1.90"* 



3.21] 
3.81"' 



9.34 
14.48 



+4.01 
2.94 
+3.63 
-3.22 
11 11+328 

26 16 +5 ' 84 
47 34+14.86 
- -13.26 
400 32 +2i2 ' 89 
1 of; ni +521.62 

J_OU.Z_L joo 1Q 



1 cc+0.52 
1 ' oo -0.34 
, 79 +0.33 
-0.23 
3+0.26 
-0.24 
3+0.32 
-0.24 
3+0.24 
-0.24 
S+0.72 
-0.73 
3 14+1-08 



1.39] 
2.10j 
1.79j 
3.68^ 



98 fiQ+28.52 
zo.u»_ 20 4 Q 
4+22.36 
1 -18.47 
-4 83+16-48 
^■°°-19.77 
31.46l 22 ;f 4 

63.93l 46 5 ; 24 

,+88.69 
1-68.50 
407 Ql+307.45 
tui -»i_276.37 



25.94"* 



53.51" 1 



3.55 
3.50 
3.52 
3.10 
3.61 
3.71 
4.09 



+0.91 
1.26 

0.8!) 
1.08 
+0.91 
1.30 
+ 1.04 
0.88 
+0.82 
1.03 
+ 0.81 
-1.25 
+0.57 
1.11 



ETG 
LTG 



23.27 
10.82 



+3.05 
-2.97 
13.05 
10.46 



1.49 
3.64 



+0.16 
0.13 
+0.82 
1.16 



41.94^ 



-10.22 
-9.29 

59.49l 42 :u l 



1.71 
3.86 



+0.29 
0.21 
+0.70 
1.04 



91 c;i+2- 85 



3.25 



+5.85 
4.51 



1.72 
2.81 



0.19 
0.16 
+ 1.36 
1.52 



19 7 o+16.79 

iy.^_ 12 41 
-26.87_i?'_5 



3.03 
3.41 



+ 1.12 
1.03 
+0.94 
1.19 



The kernel 9ggm(x) i s n °t normalised, i.e., 
J dxfem(x) 7^ 1) such that the aperture statistics as- 
sumes a value that depends not only on the underlying 
3D bispectrum B gg m but also on the normalisation. In 
order to make measurements comparable for different lens 
and source samples, we define a normalised statistics B ssm 
through the relation 

{Af 2 M ap }(9 1 -e 2 ;e 3 ) (38) 



rxh 

= Sggm (Rl,R 2 ,R 3 ) d X fam (X) 

Jo 

with Ri := /K(Xmax)#i- We emphasise that by this definition 
Bggm is not a deprojection of the angular aperture statistics 
to the spatial 3D bispectrum. This would involve the inver- 
sion of the x-integral. Instead we, in effect, normalise the 
statistic by the area J dxq gg m(x)> an d we convert angular 
scales to projected physical scales through the angular di- 
ameter distance /K(Xmax) at maximum weight g gg m(Xma.x). 
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Figure 10. Radial weight of the matter-galaxy bispectrum for the different lens samples used in this study: (_/v" 2 M ap } (top row) and 
(A/"M ap ) (bottom row). The left column is for lenses with z p hoto £ [0.2,0.44] (low-z), the right column for 2 p hoto £ [0.44,0.6] (high-z). 
L1-L6: Af r -luminositics increasing from -17.8 mag to -22.4 mag; sml-sm7: increasing stellar masses from 7 X 10 9 Mq to 4 X 10 11 Mq; 
ETG: early-type galaxies; LTG: late-type galaxies. 



The same line of reasoning can be applied to the second 
G3L aperture statistics for which we obtain 



{MM! p )(e 1 ;8 2 ;8 3 ) 

= Bmmg {Rl,R2,Rl 
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with its own radial filter 



gmmg(x) 
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Examples of kernels q , mm g (x) relevant for this work are de- 
picted in the bottom row of Fig. 10. 

By definition the aperture statistics (M n M^), as mo- 
ments of smoothed density contrasts on the sky, are dimen- 
sionless. As g ggm (x) nas the dimension [length -1 ], Eq. (37), 
we deduce from Eq. (35) that the u-filtered 6 ggm is also 
dimensionless. This becomes also obvious from Eq. (38) be- 
cause the normalisation integral is dimensionless. A similar 
argument applies to the dimensionless B mmg . 

In our analysis, we estimate the equilateral Z? ggm 
or S mmg amplitudes of all samples at a common 
comoving length scale of Ri = i?iM pc = 1 ft^pMpc (or 
k ~ \/2/-Rim P c ~ 1.4 /iiooMpc -1 for the exponential it- 
filter). For this purpose, the power-law fits in Table 3 to 
the aperture statistics are employed, which essentially de- 
scribe the data at the scales of interest, to interpolate in the 



case of S ggm 

Bggm(R) ■— B ggm {R, R, R) 



(41) 
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= :8 ggm (-BiMpc) 



We are quoting only the amplitude 6 ggm (i?iM P c) in the fol- 
lowing. Likewise for the matter-matter-galaxy bispectrum 



Bmmg we have 



Bmmg(R) '■— B 
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(42) 



x 1' 



R 



/i 10 oMpc 



/ Xh dxg mmg (x) V Kol Mpc 

— :S mmg (-RiMpc) 

To assess the uncertainty in the bispectrum amplitude, we 
marginalise over the uncertainties in Ao, the aperture statis- 
tics amplitude at 1 arcmin, and a, the power-law index, tak- 
ing into account the correlation of their errors. A value of 
fiiMpc corresponds to an aperture scale radius between 2.5 
to 3.5 arcmin depending on the mean redshift of the lens 
samples. The compiled results are plotted in Fig. 11 - one 
of the main results of our study - to highlight the trends 
with M r magnitude and stellar mass. As before, only mea- 
surements with highly significant detections are included in 
the plot. 
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Figure 11. Trends of the (significant) galaxy-matter bispectra £? m m g (-R) and £> ggm (R) for k = \/2/R = 1.03 /173MPC -1 with function of 
M r -luminosity (left panel) and stellar mass (right panel). Shown arc both redshift bins together. The left panel also includes f? ggm (ii) data 
points of the early-type galaxy sample. These data points are somewhat offset with respect to their actual mean (M r ) = — 21. 88(— 21.91) 
(high-z in brackets) for clarity. The shaded area highlights the la constraints of £? ggm (-R) for low-z and high-z combined. 



At the corresponding M r magnitude range, we include 
also the normalised (A/" 2 Af ap ) amplitude of the ETG sample. 
Their amplitude is somewhat higher in comparison to L5 and 
L6. This can be explained by the fact that the LTG sample 
is included in the L samples but not in the ETG sample of 
similar luminosity: the LTG have a normalised amplitude 
considerably smaller than that of the ETG. 

5.4 Third-order galaxy biasing 

SW05 introduced a set of third-order galaxy biasing param- 
eters &3,ri,r 2 to parametrise the galaxy- matter bispectra 
relative to the matter bispectrum S mmm , 



B ggm (fei,fe 2 ,x) 
-B mmg (fci,fc 2 ,x) 



&3 r 2 S mmm (fei, fe 2 ,x) 
h^r\ B mmm (fei, fc 2 ,x) 



(43) 
(44) 



The coefficients &3,ri,r 2 are also functions of fci,fc 2 ,x, 
which has been omitted here to save space. For galax- 
ies faithfully tracing the underlying matter distribution 
one finds n = r 2 = 63 = 1 for all scales. This parametrisa- 
tion generalises the earlier similar second-order galaxy bias 
parametrisation (e.g., Pen et al. 2003; Hoekstra et al. 2002; 
Tegmark & Bromley 1999) to the third-order. Our nor- 
malised G3L measurements of S ggm (-R) and Bmmg(-R) can 
be utilised to constrain the ratio r\j{ribz) for a physical 
scale R by considering the combined statistics 

^5mm g (-^) 



(45) 



Bggm(^) 

Jo Xh dxq ggm (x) (AfQ (A/k 1 ^)) 

Jo Xh dX<? mmg (x) (A/" 2 M ap ) (fl/K^Xmax)) 

for a lens sample. We assume here that both kernels q ggm 
and g mmg peak at the same distance Xmax, which is approxi- 
mately valid for our study. This bias parameters in ri/(r 2 &3) 



are smoothed in k and x with maximum weight at k « y2/R 
(equilateral triangles) and Xmax- The exact smoothing ker- 
nels are given in Appendix D. A deviation of ty(R) from 
unity hence indicates a biased galaxy population. 

We calculate the ^(R) statistics for the samples L1-L6, 
samples sml-sm5, and the ETG sample (all low-z only) for 
angular scales between one and ten arcmin. The remaining 
measurements are too noisy for useful constraints. Fig. 12 
summarises these novel measurements. The error distribu- 
tions of the ratios ^(R) are estimated by employing Monte- 
Carlo realisations of S ggm and B mmg (assumed Gaussian); 
depicted are the mean and variances 0* in the resulting 
distributions. Alternatively, one could utilise the analytic 
probability distribution function given in Hinkley (1969). 



6 DISCUSSION 

We performed a G3L analysis of approximately 100 square 
degree of the CFHTLenS data set. The data is endowed 
with photometric redshifts of galaxies and lens&t estimates 
of the PSF-corrected source ellipticities. For the first time, 
the signal-to-noise of the lensing data is sufficient to mea- 
sure third-order galaxy-galaxy lensing as a function of lens 
luminosity, stellar mass and galaxy type. The work of Simon 
et al. (2008), analysing the RCS1 data, demonstrated that 
G3L measurements are principally possible with contempo- 
rary lensing surveys. This is confirmed by this study. We fur- 
ther subdivided the lens samples in M,-luminosities, stellar 
masses, SED types, and a "low-z" (0.2 ^ z p hoto < 0.44) and 
a "high-z" (0.44 ^ z p hoto < 0.6) redshift bin by utilising the 
photometric redshifts of the lenses. We presented the G3L 
measurements in terms of aperture statistics that probes the 
angular bispectrum of the (projected) matter-galaxy three- 
point correlations. In one case ("lens-lens-shear"), the mea- 
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Figure 12. Results of the ^-statistics which probes the third-order galaxy biasing parameters of our lens samples. Plotted is ri/{b^T2) 
as function of scale k = \/2/R = \/2/(/K(Xmax)9) for the angular range 1' ^ 6 < 10'. Used are only significant measurements in the 
low-z samples. Luminosity samples are in the left panel, stellar mass samples and the ETG sample in the right panel. Data points with 
uncertainties larger than > 0.8 are not shown. The shaded area indicates the mean and standard deviation of the mean of combined 
samples. L1-L6: M r -luminosities increasing from -17.8 mag to -22.4 mag; sml-sm5: increasing stellar masses from 7x 10 9 Mq to 10 11 M ; 
ETG: early-type galaxies. 



surements quantify correlations between two lens positions 
and the surface matter density around the lens pair; this can 
be interpreted as excess surface mass density about galaxy 
pairs (Simon et al. 2012). In the other case ("lens-shear- 
shear"), it expresses correlations between a lens position 
and the surface mass density in two different directions close 
to the lens. The here adopted G3L aperture statistics has 
the practical advantage to separate E- and B/P-modes from 
these measurements, which is utilised to detect signatures of 
possible systematics in the data. On this level, no significant 
G3L systematics signals were detected. 

To reduce the impact due to intrinsic alignments of 
sources, we separated lens and source galaxy samples phys- 
ically from each other by exploiting the photometric red- 
shifts in the survey. We showed that in the ideal case of 
no radial overlap, neither Il-correlations nor Gl-correlations 
contribute to the correlator. Owing to the uncertainty in the 
galaxy redshifts, however, perfectly non-overlapping distri- 
butions are hard to achieve. We found that our low-z lens 
samples have still a small overlap of ~ 4% with the source 
sample, the high-z samples a moderate overlap of ~ 12% 
(overlap of redshift probability distribution functions). Be- 
cause of the small overlap, at least for the low-z samples we 
do not expect significant contributions from intrinsic align- 
ment correlations. To test the degree of contamination by 
GI- and Il-correlations, we compared the aperture statistics 
of the combined sm3-sm5 samples, both low-z (zd = 0.35) 



and high-z (z^ = 0.51), for two cases. In one case, we se- 
lected sources by 0.65 ^ z p hoto < 1.2 (z s = 0.93) as before. 
In the second case, we were more conservative by selecting 
only sources within 0.8 z p hoto < 1-2 (z s = 1-02), thereby 
discarding about one third of our sources. However, the lat- 
ter case reduced the overlap from 3.3% (10.9%) to 0.6% 
(2.5%) for the low-z (high-z) sample. The statistics were 
normalised by J dx?g g m(x) (/ dx<Zmmg(x)) to compensate 
the signal change owing to the different source depths. The 
maximum signal increase is ~ 30% for (JVM^ p ) of the high- 
z lens sample. Fig. 13 shows the difference in normalised 
statistics for fixed lens samples but varying source depths. 
Here we assumed that /K(Xmax) is identical for both com- 
pared signals, i.e., both signals were differenced at the same 
aperture scale radius. This assumption is accurate within a 
few percent here. As expected, the difference is consistent 
with zero, the level of GI/II-systematics in the statistics is 
therefore negligible within the measurement errors. 

As second possible source of systematics we identified 
the magnification of the lens number densities by matter 
fluctuations in front of the lenses. To first order this effect 
is comparable to the aperture mass skewness (Mf p (#)} asso- 
ciated with sources at redshift z s ~ 0.4. We estimated this 
effect to be of the order of<10~ 8 at#~l arcmin. We 
therefore conclude in comparison with our measurements 
that in the range V < 8 < 10' the lens number density 
magnification effect is negligible for both (A/" 2 M ap ) and for 
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Figure 13. Plotted are the differences of 
(A/" 2 M ap )(0)/ f dxijggm(x) f° r sources with z photo > 0.65 and 
^photo > 0.8, and the analogue for (A/"Mf p )(0)/ / <?mmg(x)- 
For fixed lens samples (sm3-5 combined; either low-z or high-z), 
the result is expected to be consistent with zero. Note that 
neighbouring error bars are correlated. 

(JVMfp). For a more detailed investigation of systematics, 
however, realistic models of intrinsic alignments and magni- 
fication for third-order correlations are required, which are 
currently unavailable. 

The observed aperture statistics depend on the redshift 
distributions of lenses and sources. Results of the statistics 
would hence differ when changing the source or lens red- 
shift distribution, even if the underlying comoving spatial 
3D matter-galaxy bispectrum were unchanged throughout 
the light-cone. In order to partially correct for this effect, 
we employed a new technique that normalises the aper- 
ture statistics with the lensing efficiency and relates angu- 
lar scales to effective spatial scales; this yields two kinds of 
galaxy-matter bispectra B gsm (R) and Z3 mmg (7?), originating 
from (AT 2 M a p} and (JVM 2 p ), respectively. The normalised 
statistics obtained are basically band bispectra due to the 
smoothing of the exponential ti-filter in fc-space and the ra- 
dial smoothing of the lensing kernels. Only by means of our 
or similar normalisation schemes, measurements for different 
lens samples or same galaxy selections at different redshifts 
become comparable. In particular, the application of our 
normalisation simplifies the comparison with results from 
future studies. The problem of unnormalised measurements 
becomes particularly obvious for {AT 2 M ap } of the low-2 LI 
sample in Fig. 6 (top left panel) in comparison to the L6 sam- 
ple in the same panel: for 9 > 2 arcmin both measurements 
are basically identical, although the normalisation reveals 
that the lower luminosity galaxies have a smaller bispec- 
trum amplitude (Fig. 11, left panel). This effect results from 
a completely different redshift distribution of the LI lenses 
that, due to sample incompleteness, have a mean redshift of 
z ~ 0.22 instead of L6's z ~ 0.34. 

We estimated the measurement errors directly from the 
data by Jackknifing the signal variance across the survey 
pointings. Ideally, with statistically independent pointings 
this would properly account for uncertainties due to source 
shape noise, sampling noise and cosmic variance. However, 



Table 4. Values of the normalised galaxy-matter bispectra 
-Bggm(K) (top half) and B mmg (ii) (bottom half) at R = 
1 /i^qqMpc for the late-type (LTG) and early- type (ETG) galaxy 
sample. 

Sample low-z high-z 



ETG 
LTG 


+l.04±°;i? 

+ 4.37±^ 


x 10" 
x 10" 


4 

6 


+9.48:! 
+1.07j 


-1.19 
-1.23 
-1.67 
-1.02 


X 10" 
X 10" 


-5 
-5 


ETG 
LTG 


+6.93i^ 


x 10" 
x 10" 


-5 
-6 


+1.42 j 
-1.39j 


-1.37 
0.92 
-0.95 
-1.62 


x 10" 
x 10" 


5 
-5 



the pointings are bundled together in large fields W1-W4 
with extensions of several degrees across the sky. This makes 
pointings of the same field partly correlated. Therefore, the 
cited uncertainties are probably somewhat too optimistic in 
the sense that they underestimate the cosmic variance. 

To refine the previous RCS1 measurement in Simon 
et al. (2008) for different galaxy populations and to investi- 
gate the dependences of bispectra amplitudes on galaxy pop- 
ulations, we focussed here on equally-sized apertures. This 
gives most weight to the equilateral bispectra. We found 
that the aperture statistics are reasonably well described by 
a power law on angular scales ranging from roughly two ar- 
cmin to ten arcmin. On smaller angular scales, we observe 
evidence for a change of slope, but we are also increasingly 
affected by the transformation bias. For instance, (A/" 2 M ap ) 
of the fainter low-z L4 sample clearly flattens below ~ 2 
arcmin. Qualitatively, this behaviour is also observed in the 
semi-analytic galaxy models studied in Saghiha et al. (2012), 
see their Fig. 8. A similar change of slope, maybe also a 
steepening, is visible for the {AfM 2 p ) statistics of the lu- 
minosity samples. A comparison of galaxy models to our 
measurements requires a careful replication of the galaxy 
sample selections, their uncertainties, and the survey incom- 
pleteness. Moreover, as concluded in Saghiha et al. (2012), 
no reliable galaxy model is currently available to predict 
the correct amplitude of G3L measurements - or to even 
double-check whether our results may be strongly effected 
by galaxy selection effects. We hence postpone this task to 
a future paper. 

The measurements of S ggm utilising ETG and LTG 
pairs (subdivision of the combined L5/L6 sample) show that 
the excess mass around pairs is a strong function of galaxy 
type. The high excess mass signal of ETG is comparable 
to the strong signal of pairs in our sm-samples with stel- 
lar masses of ~ 10 11 Mq, whereas the excess mass of LTG 
is consistent with zero in our case (Table 4). A plausible 
explanation for this is the fact that many early-type galax- 
ies in the ETG sample ((M*) ~ 6 x 10 10 M©) are satellites 
in dense cluster environments, whereas LTG are frequently 
field galaxies. This was, for instance, found by the GGL 
study of Mandelbaum et al. (2006) in SDSS. The splitting 
into ETG from over-dense and under-dense regions that was 
conducted in this study is actually comparable to the lens- 
lens-shear correlation function because the G3L correlator 
gives more weight to pairs in cluster environments, sim- 
ply because more pairs are found in these regions. Recently, 
Saghiha et al. (2012) studied the excess mass for two state- 
of-the-art semi-analytic galaxy models. Although the G3L 
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amplitudes and colour distributions of the two considered 
models are inconsistent, both models predict a large differ- 
ence in (M 2 Map) for za = 0.17 red and blue galaxies up to 
a factor of ~ 10 2 at 8 ~ 2 arcmin. With our uncertainties in 
the LTG signal, we estimate the difference to be at least a 
factor ~ 10, strongly confirming the previous prediction. 

By forming ratios of normalised bispectra, our new 
statistic ^(R) := B m mg(R) / 6 ggm (i?) approximately yields 
the ratio 7-1/(^63) of (smoothed) third-order biasing param- 
eters (SW05). The details of the smoothing are determined 
by the shapes of the it-filter, the peaked kernels q mmg , g ggm , 
and to some extent also the matter bispectrum (Appendix 
D). Deviations of from unity indicate galaxies that not 
faithfully trace the underlying matter density field, i.e., bi- 
ased galaxies. This new technique for investigating galaxy 
bias with lensing advances the methodology of van Waer- 
beke (1998) and Schneider (1998) that focus on second-order 
galaxy bias. The application of the latter found that galax- 
ies are generally biased tracers (Hoekstra et al. 2001; Simon 
et al. 2007; Julio et al. 2012). We confirm this finding by 
employing third-order galaxy-matter correlations. 



7 CONCLUSIONS 

• We detect G3L with unprecedented high significance in 
the CFHTLenS for galaxy populations of different luminos- 
ity, stellar mass, and SED type. This applies to both third- 
order galaxy-galaxy- matter correlations (6 ggm ) and galaxy- 
matter- matter correlations (Z? mmg ). 

• We find that the (equilateral) galaxy-matter bispec- 
tra are, within the remaining statistical errors, reason- 
ably well scale-invariant for the spatial (comoving) scales 
0.3 Mpc -1 < k < 2.2 Mpc -1 . On smaller scales, not included 
in our power-law fits, there are indications of deviations from 
the power-law shape in several cases. 

• The low-z and high-z counterparts of the same lens 
samples yield very similar bispectra amplitudes S ggm (Fig. 
11) and slopes (Fig. 9). This points to little evolution of the 
bispectrum between redshift z ~ 0.3 — 0.5, especially for our 
Mr-selected galaxies. There is, however, some evidence for 
a change in the amplitude of Z? ggm for stellar-mass selected 
galaxies below ~ 10 11 M©: high- 2 lenses show a lower am- 
plitude than the low-z lenses (right panel of Fig. 11). This 
implies an increase of excess mass about pairs of galaxies 
of fixed stellar mass with time, as, e.g., may be expected in 
a CDM scenario due to the continuous accretion of matter. 
Evolution trends of Z3 mmg are unclear due to the measure- 
ment uncertainties in the high- z samples. 

• For Z3 ggm the slope and the amplitude is a changing 
function of galaxy luminosity, stellar mass and galaxy type. 
The amplitude change is also observed for Z3 m m g . Brighter 
or more massive galaxies (by stellar mass) exhibit a steeper 
bispectrum, which implies that the excess mass is more con- 
centrated in these cases. Moreover, there is a clear trend to- 
wards higher amplitudes for both more luminous and more 
massive galaxies. This shows that more luminous or massive 
galaxies, or galaxy pairs in the case of S ggm , inhibit denser 
environments than fainter or lighter galaxies. 

• We observe a strong signal for the excess mass around 
early-type galaxies (ETG) pairs. Late-type galaxy pairs 
(LTG), on the other hand, have a signal that is consistent 



with zero when studied as aperture statistics. This remark- 
able observation is in excellent agreement with the recent 
prediction of Saghiha et al. (2012) based on semi-analytic 
galaxy models. The measurement therefore suggests that 
virtually all signal in this magnitude range originates from 
ETG pairs, and possibly mixed pairs of ETG and LTG, 
rather than from LTG pairs. This can be explained by the 
fact that a large fraction of ETG are satellite galaxies in 
cluster. By explicitly mapping out the excess mass around 
LTG and ETG galaxy pairs we have also found that both 
maps are fundamentally different in their amplitudes as well 
as in their general appearance. 

• The mismatch between S ggm and Z? mmg for the same 
lens galaxy sample immediately indicates galaxies biased 
with respect to the matter distribution. This mismatch is 
captured by the galaxy bias statistics ^(R) (Fig. 12) that 
shows for our low-z samples values comparable for a wide 
range of galaxy luminosities and stellar masses. Therein, 
we probe the non-linear regime on scales smaller than k ~ 
0.8 Mpc~ . We find best constraints on ^(R) with the stellar 
mass samples, which has for all samples sml-sm5 and scales 
combined (minimum- variance weighted) an average value of 
\t = 0.51 ± 0.07. This shows - for the first time employ- 
ing third-order lensing statistics - that galaxies are biased 
tracer of the matter density field. Although ^(R) indicates 
that the ratio ri/(r2&3) stays relatively constant with scale, 
with a possible shallow local minimum at k « 1.8 Mpc -1 , the 
additionally observed change of the bispectrum amplitudes 
with galaxy luminosity or mass means that the individual 
bias parameters have to differ between the galaxy samples. 

• Finally, we emphasise that theory is lacking behind in 
interpreting the G3L measurements. Reliable model predic- 
tions , e.g., in the vein of Takada & Jain (2003), are needed, 
not only to properly interpret the measurements, but also to 
gain a better understanding of systematics and to verify that 
selection effects in the data do not spoil the measurement. 
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APPENDIX A: MULTIPLICATIVE SHEAR 
BIAS 

Miller et al. (2012) discusses a calibration scheme for corre- 
lation function estimators involving shear estimates from the 
lens&t pipeline. For details, we refer the reader to the men- 
tioned article, Sect. 8.3 and 8.4. Analogous to the calibration 
scheme of the two-point shear-shear correlation function de- 
tailed therein, we divide G cst , Eq. (16), and G±\ Eq. (17), 
by l + A'g;($i,i92,</>3) and 1 + K G± (#i, #2, ^3), respectively. 
Both calibration factors are given by 

1 + Kg(#i,& 2 ,cf>3) (Al) 

N d N d N s 

E E E^a + ^K 1 * 203 

_ i = l j = l fc = l 

_ N d N d N s ' 

E E E ^K 1 * 3 
i=i 3=1 fc=i 

1 + Jf 0± (0i,0a, 0a) 
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N A JV B N s 

E E E ^^(l + m,)(l + m,)A^ 
t=l j=l k=l 

N d JV 3 JV B 

E E E wA***' 



i=l 3=1 k=l 



ijk 



The multiplicative bias factors rru, provided in the 
CFHTLenS catalogue 3 for each source, are functions of the 
source signal-to-noise and angular size. Note that for non- 
vanishing values of 1 + mm > 0, the calibration is mathe- 
matically equivalent to employing the transformation i— ► 
6j/(l + T7ij) and u>i H> + m,) in the estimators C/ est and 



APPENDIX B: LENS SAMPLES SUPPLEMENT 
Bl Angular clustering of lenses 

The angular correlation function lj(8) of the lenses as a func- 
tion of separation 9 is approximated by a simple power law 
(Peebles 1980): 



oj(8) = A* 



IC 



(Bl) 



where A u is the clustering strength at a separation of one ar- 
cmin, A the slope of the correlation function and the constant 
offset IC, the integral constraint (Groth & Peebles 1977). We 
find this fitting function to be a good description of ui(0) be- 
tween 0'.2<9< 10'. The estimator of w(0) by Landy & Sza- 
lay (1993), employed for this paper, requires the repeated 
counting of galaxy pairs in separation bins for random real- 
isations of unclustered galaxy distributions, factoring in the 
incompleteness of the survey. The figures quoted in Table 1 
consider the masks of individual survey pointings, but pre- 
suming the same survey completeness within regions that 
are not masked out. For the final to(0), all pair counts from 
all individual survey fields are added so that, in effect, fields 
with more galaxies attain a higher weight in the average. 
The binned oj(0) is stored as data vector d. 

Pointing-to-pointing Jackknife sampling To esti- 
mate the statistical uncertainty of d, we prepare a set of N p 
Jackknife samples df 1 , where di is the combined data vec- 
tor omitting the ith patch. The Jackknife covariance of the 
sample mean is then (Shao 1986; Norberg et al. 2009): 



i p t 



where 
Adf:: 



(7Vp-l)(d-dn 



(B2) 



(B3) 



For Table 1, a power law fit is applied to the ensemble av- 
erage d of all pointings, taking into account the Jackknife 
covariance C. The ^-binning in d is also applied to the power 
law model, Eq. (Bl), by averaging the model over the width 
of each bin. Note that the inverse of C, required for the 
likelihood analysis of the model fit parameters, has to be 
corrected to obtain an unbiased estimator of the inverse co- 
variance (Hartlap et al. 2007). Similar corrections of inverse 

3 Publicly available under http : //www . cadc-ccda.hia-iha. 
nrc-enre .gc . ca/community/CFHTLens /query .html 



covariances are applied throughout the paper without fur- 
ther mentioning. 



B2 Completeness of lens sample 

Here we define a parameter that quantifies the completeness 
of our lens samples. First we define the comoving volume 
V(z\,Z2) confined by the redshift boundaries z\ ^ z < Z2, 



v( Zl , Z2 ) = n 



X(zi) 



xOo) 



dx/ K (x) , 



(B4) 



where fl is the solid angle of the patch field-of-view minus 
the solid angle of mask regions, and 



dz' 
E(z>) 



(B5) 



where E(z) := H(z)/Ho is the Hubble parameter H(z) as a 
function of redshift normalised to Hq. 

Due to the incompleteness in our flux-limited survey, a 
galaxy is only visible up to a certain redshift z ma x- In gen- 
eral, and especially for faint galaxies, one can expect the 
limit 2 m ax to be a complicated function of intrinsic galaxy 
properties and position, survey instrumentation, survey con- 
ditions and the data reduction pipeline. Nevertheless, here 
we take the simplistic view that the main factor is the ap- 
parent i'-band magnitude of the lens (extinction corrected) , 
which is limited to il ^ 22.5, such that our lens samples are 
predominantly magnitude limited. We further assume that 
a if-correction is negligible over the redshift bin \z\ , 22] of 
interest. Under these circumstances, one finds implicitly for 



D L (z maX ) = 10' 



0.4(m H , 



-m) 



D^z) 



(B6) 



where z is the redshift of the galaxy, m its i'-band magni- 
tude, and mumit = 22.5 the asserted magnitude limit of the 
lens catalogue. By Dl(z) = (1 + z)fK(x( z )) we denote the 
comoving luminosity distance as a function of redshift. 

In order to quantify for Table 1 the completeness of 
a sample of N s galaxies, we estimate over which fraction 
/ c = V(zi, z max )/V(zi, 22) an observed galaxy in the sam- 
ple would be observable. We take the average of all volume 
fractions of all lenses in a sample, 



1 - \ / 2 



i=l Jz l 



V^z^MINOzw,*;^)) 



V(z u za) 



(B7) 



and marginalise over the uncertainties in the galaxy red- 
shifts, quantified by the p.d.f. Pi(z)dz. Importantly, 2 maXii 
denotes the maximum redshift at which the ith galaxy would 
still be included within the galaxy catalogue, complying with 
all survey and sample selection criteria. A completeness pa- 
rameter close to unity means that essentially all galaxies in 
the sample are visible throughout the entire volume, whereas 
f c <C 1 indicates a significant fraction of galaxies that is only 
visible in a small subvolume at lower redshift. Obviously, / c 
is merely an estimator (upper limit) for the sample com- 
pleteness as galaxies already not observed at redshift z\ are 
not accounted for. Note that the solid angle SI cancels inside 
the expression for / c and is hence not needed. 
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Figure CI. Plots of (A/"A/ x M ap ), {MM"]) (LI low-z and sm7 
high-z) in comparison to the E-mode (lines). Both samples failed 
the 95% confidence level null test for the P/B-modcs, sec Table 
2. 



APPENDIX C: SYSTEMATICS 



C2 Field dependence of the G3L signal 

The CFHTLS wide survey consists of four contiguous fields 
Wl (~ 72deg 2 ), W2 (~ 33deg 2 ), W3 (~ 49deg 2 ) and W4 
(~ 25deg 2 ); the field areas do not include masking or ex- 
cluded fields due to significant PSF residuals. The fields are 
well separated on the sky and were observed at different 
times of the year. By splitting the measurements into WT 
W4, we check whether the G3L measurements are compa- 
rable for different subsets of the data. To have a possibly 
large sample for this test, we combine the signals of the 
samples LI to L6 for each field, see Fig. C2. Only mea- 
surements from pointings within the same fields W1-W4 are 
combined, their statistical uncertainties originate from the 
Jackknife technique (as in Eq. B2). Therefore, the error bars 
do not include the cosmic variance between the fields, which 
should be most prominent at the larger angular scales. 

We find excellent agreement between the measurements, 
considering that statistical uncertainties at larger scales are 
higher than indicated and that errors between neighbouring 
angular bins are correlated. In particular, this separation of 
data shows that the G3L signal does not originate from sin- 
gular fields that are extreme outliers compared to the others. 
Since the uncertainties of the final combined measurements 
are based on the pointing-to-pointing variance of the entire 
survey, a possible systematic deviation of one field will be 
included as systematic error inside the error bars. 



CI B/P-mode consistency with null 

As indicator of possible systematics in the estimators, 
we test (Af 2 M ± ) and the combined (AfM ± M ap ), (A/"M|> 
against the null hypothesis. A null measurement should re- 
sult in a 

A x 2 =c?C 1 d, (CI) 

that is statistically consistent with a vanishing signal, with 
d being a vector consisting of the measurements for the P- 
mode ((A/" 2 Afa P )) or both the P- and B-mode ((A/"M 2 p )). By 
C we denote the Jackknife covariance of the measurements 
as obtained from the variance of B/P-mode measurements 
in the pointings, as explained in Sect. Bl. This covariance 
is larger than a null hypothesis covariance as it possibly 
also contains power from B/P-modes present in the data. A 
true null model would contain only power from galaxy shape 
noise and sampling noise. The test results can be found in 
Table 2. Measurements inconsistent with a null signal (95% 
confidence) are underlined, thus for Ax 2 ^ 2.0 (1.68) per 
degree of freedom for (A/" 2 A/ ap > ((A/"Mf p )). In total, we find 
two lens samples that fail the test; they are plotted in com- 
parison to their E-mode in Fig. CI. In both cases the failures 
are related to the (A/"M 2 p ) statistics and significantly neg- 
ative B-modes. Note that errors between neighbouring bins 
are strongly correlated. 

Finding two measurements out of in total 57 that fail 
the 95% test is what we would expect as false positive rate. 
We therefore conclude that the influence of systematics on 
the E-mode measurement that reveal themselves via the P- 
or B-modes is likely to be small compared to our measure- 
ment uncertainties. Note that the sm7 sample that failed 
the null test is not used in the final analysis because the 
corresponding E-mode signal is consistent with zero. 



APPENDIX D: THIRD-ORDER GALAXY 
BIASING 

The values 6 ggm (S m mg) measure the it-filtered bispectrum 
B gS m (5 m mg), radially smoothed with maximum weight at 
Xmax- The maximum weight of the u-filter in Fourier space 
is at k = y/2/R for a given real space scale R. From the 
definition of the ^/-statistics, Eq. (45), from Eq. (34) and 
from a similar equation for (NM 2 p ) it follows that 

(Dl) 

_ / dxd 2 fcid 2 fc 2 F mmg (ki, k 2 , x; R) (b 3 ri) (fei, fca,x) 
/dxd 2 fcid 2 feF ggm (fci,fc2,x;#) (blr 2 )(ki,k 2 ,x) 
where the smoothing kernels in (fe, x)-space are 

F jnmg (ki,k 2 ,x;R) (D2) 

= gmmgfa) x 

fo h dXgmmg(x) 

u(kiA)u(k 2 A)u(\ki + fe2|A)B mmm (fei, fe2,x) > 
F egm (k 1 ,k 2 ,x;R) 

= g gg m(x) x 

u(fciA)u(fc 2 A)w(|fei + fe 2 |A)_B mmm (fc 1 ,fc 2 ,x) 

with A := -R/k(x)//K(Xmax). As can be seen, the detailed 
weight within a band defined by the width of the u-filter is 
also determined by the actual matter bispectrum -B mmm . 

We can further exploit the statistical isotropy of the 
galaxy-matter bispectra, which means that both -Bmmm and 
the galaxy biasing parameters ri,r 2 ,bz are only functions 
of |fei|, \k 2 \, cj>; (j> i s the angle spanned by fci and k 2 . The 
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Figure C2. Measurements of the aperture statistics (top row: (A^ 2 M ap }, bottom row: (J\fM^ p )) for the combined L1-L6 sample (low-z 
and high-z rcdshift bin separately). The measurements are split for the four fields W1-W4. The error bars indicate only the pointing-to- 
pointing variance within the same field. The fields vary in size and thus the sizes of lens and source catalogues vary. Lines connect the 
data points to guide the eye. 

previous expressions therefore simplify to 

9(R) (D3) 
_ / dxd(t>dk 1 dk 2 k 1 k2 F mmg (- . .) (6 3 ri) (fci, ft 2 , 0, x) 
/ dxd(/>dk 1 dk 2 k 1 k2 F gsm (. . .) (6§r 2 ) (fci, k 2 , <j>, x) 
where 

F mmg (...) (D4) 

= gmmg(x) x 

u(k 1 A)u(k 2 A)u(\ki + fc 2 |A)B mmm (fc 1 , k 2 , (f>, x) , 
F ggm (. . .) 

_ 9ggm (x) 

"~ / Xh dXfeg m (x) 
u(kiA)u(k 2 A)u(\k 1 + fc 2 |A)B mmm (fci, k 2 , <f>, x) , 

and 

|fei + k 2 \ = \Jkf + kl + 2fcifc 2 cos0 . (D5) 

Note that for equilateral triangles we have ki — k 2 = |fci + 
k 2 \ and thus cos0 = —1/2. 



